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BLAZARS AS ULTRA-HIGH-ENERGY COSMIC-RAY SOURCES: IMPLICATIONS FOR TEV GAMMA-RAY 
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ABSTRACT 

The spectra of BL Lac objects and Fanaroff-Riley I radio galaxies are commonly explained by the one-zone 
leptonic synchrotron self-Compton (SSC) model. Spectral modeling of correlated multiwavelength data gives 
the comoving magnetic field strength, the bulk outflow Lorentz factor and the emission region size. Assuming 
the validity of the SSC model, the Hillas condition shows that only in rare cases can such sources accelerate 
protons to much above lO' eV, so > 10 eV ultra-high-energy cosmic rays (UHECRs) are likely to be heavy 
ions if powered by this type of radio-loud active galactic nuclei (AGN). Survival of nuclei is shown to be 
possible in TeV BL Lacs and misaligned counterparts with weak photohadronic emissions. Another signature 
of hadronic production is intergalactic UHECR-induced cascade emission, which is an alternative explanation 
of the TeV spectra ofsome extreme non-variable blazars such as lES 0229H-200or lES 1101-232. Westudythis 
kind of cascade signal, taking into account effects of the structured extragalactic magnetic fields in which the 
sources should be embedded. We demonstrate the importance of cosmic -ray deflections on the 7-ray flux, and 
show that required absolute cosmic-ray luminosities are larger than the average UHECR luminosity inferred 
from UHECR observations and can even be comparable to the Eddington luminosity of supermassive black 
holes. Future TeV 7-ray observations using the Cherenkov Telescope Array and the High Altitude Water 
Cherenkov detector array can test for UHECR acceleration by observing > 25 TeV photons from relatively 
low-redshift sources such as lES 0229H-200, and > TeV photons from more distant radio-loud AGN. 
Subject headings: galaxies: active — gamma rays: galaxies — cosmic rays 



1. INTRODUCTION 

Active galactic nuclei (AGN) with extended radio jets pow- 
ered by super-massive black holes are among the most lumi- 
nous objects in the low-redshift universe. Since 2004, when 
the present generation of imaging atomspheric Cherenkov 
telescopes began to operate, the number of extragalactic 
sources detected at > 0.1 TeV (very -high energy; VHE) en- 
ergies has grown rapidly, and is nearly fifty.^ The Fermi 
Gamma-ray Space Telescope, now in its fourth mission year, 
is providing a wealth of new discoveries on 7-ray galax- 
ies. In the high-confidence clean sample of active galac- 
tic nucl ei associations in t he First Fermi LAT AGN catalog 
(lLAC: lAbdoetal. 2010ah . more than 600 7-ray blazars, di- 
vided about equally into BL Lac objects and flat spectrum 
radio quasars (FSRQs), were reported. New classes of GeV 
7-r ay galaxies, e.g., r adio-loud narrow line Seyfert galax- 
ies dAbdo et al. 2009al) and star-formin g galaxies p owered by 
supernovae rather than black holes dAbdo etal. 2010b,) . fol- 
lowing clo sely the VHE det ections of the starburst galaxie s 
NGC 253 flAcero et al . 2009*) and M82 ( lAcciari et al. 2009cl) . 
are now firmly established. Moreover, the spectral energy dis- 
tributions (SEDs) of radio galaxies detected at GeV and VHE, 
being misaligned by large (> 10°) angles to the jet axis and 
thought to be the par ent population of blazar s in geometrical 
unification scenarios dUrrv & Padovani 1995h . are helping to 
reveal the blazar jet g eometry. About 10 such sources are now 
detected with Fermi (lAbdo et al. 2010cl) . 
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Radio-loud AGN detected at TeV energies consist 
mainly of high-synchrotron-peaked BL Lac objects, in- 
cluding the ultra-v ariable TeV blazars Mrk 421 (z=0.031; 
iFossati et al. 2()08h . Mrk 501 (z=0.033; Albert et al. 20071^ 
and PKS 2155-304 (z=0.116; .Aharonian et al. 2007al) . 
and the appa rently non-variable TeV blazars lES 
0229- H200 (Ahai-onianetal. 2007b) and lES 1101- 
232 d Aharon ian et al. 2007c). Extragalactic VHE 7-ray 
galaxies include several Fanaroff-R iley (FR) class I radio 
galaxies (Cen A, M87 , NGC 1275;rAharonianetal. 2009at 
lAharonian et al. 2006t lAleksic et al. 201 Id respectively 



Cen A dAbdo et al . 2010d'),^187 (Abdo et al. 2009c 
and NGC 1275 (JAb do et al. 2009bl; Kataoka et al. 20 1 (| 
iBrown & Adams 201 ll) have also been detected at GeV 
energies. This hst also include s the head-tail radi o 
galaxy IC 310 dAleksic etal. 20101: iNeronovetal. 2010h . 



intermediate-synchrotron-peaked 
3C 66A (A cciari et al. 2009i 
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Aliuetal. 2009: 



Abdo etal. 2011 



Abdo etal. 2011 



BL Lac dAlbert etal. 2007a : 

the Ge V luminous, high-redshift 

FSRQs 3C 27 9 (z=0.538; Aleksic' etar2011ah . PKS 1510- 
089 (z=0.361; 'Wagner etal. 2010), and 4C H-21.35 (PKS 
1222h-216,z=0.432; Aleksic et al. 2011b). 

The 7-ray data from weak-lined BL Lac objects and 
FR-I radio galaxies are generally well fit with the 
standard nonthermal electronic synchrotron self-Compton 
(SSC) relativistic jet model (e.g..lMastichiadis & Kirk 19971: 
iTavecchio et al. 19981: iKato et al. 20061) . but the use of 
archival data for highly variable blazars gave large parameter 
uncertainties in the past. With the recent simultaneous multi- 
wavelength data sets for many sources, accurate parameter es- 
timation can be made, either from simple scaling results in the 
Thomson regime, or from detailed spectral calculations taking 
into account the Klein-Nishina (KN) effect that is relevant for 
high, UV/X-ray synchrotron-peaked blazars. 

The paper is organized as follows. In Section 2, we assem- 
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ble the derived parameter values obtained in various analy- 
ses of typical BL Lac objects and FR-I radio galaxies. From 
these numbers, we obtain maximum energies of cosmic rays, 
and show that protons can be accelerated to > 10 EeV en- 
ergies only in a few radio galaxies and flares of BL Lac ob- 
jects. Then we explore the associated hadronic signatures ex- 
pected from TeV blazars in the case where jets of BL Lac 
objects and FR-I radio galaxies are accelerators of ultra-high- 
energy cosmic rays (UHECRs, with energies above the an- 
kle of w 10 ■ eV), assuming the validity of the leptonic 
SSC parameters inferred from rapidly variable AGN. We dis- 
cuss observable signals produced in the source and calculate 
those generated outside the source from both 7 rays and UHE- 
CRs escaping from the jet accelerator and passing through the 
~ Mpc-scale regions of cosmic structure, magnetized clus- 
ters and filaments, and the larger ^100 Mpc-scale voids of 
intergalactic space. For some apparently non-variable TeV 
blazars, where the one-zone synchrotron/SSC model typi- 
cally requires extreme parameters for fits, the cascade radi- 
ation can be a crucial component of the high-energy radiation 
spectrum, depending on the strength of the intergalactic mag- 
netic field (IGMF) in voids, as has recently been proposed 
to explain the 7- ray spectra of some extreme blazars such 
as lES 0229-H200 (lEssev & Kusenko 2010l:lEssey et al. 20101: 
lEssevet al. 2011 *). We focus on such cascade emissions in 
the VHE range in Section 3. Using numerical calculations, 
we also demonstrate the importance of structured extragalac- 
tic magnetic fields (EGMFs) in clusters and filaments for fu- 
ture 7-ray detectability by the Cherenkov Telescope Array 
(CTA) and the High Altitude Water Cherenkov (HAWC) de- 
tector array. Implications of this study and a summary are 
given in Section 4. Throughout this work, the cosmological 
parameters are taken as //q = 71 km s"' Mpc~', ^^ = 0.3, and 
r2A = 0.7. 



2. VHE BLAZARS AND UHECRS 

Of the wide variety of source classes that could po- 
tentially accelerate_UHECR5^_^iKluding^_^oi__ex^^ 
GRBs (IWaxman 19951; IVietri 19951; iMurase et al. 2008al) . 
fast rotating magnetars (lArons 2003h, structure forma- 
tion shocks in galaxy clusters JNorman et al. 19951 ; 
iKang et al. 1996t llnoue et al. 2007h, and quasar rem- 
nants ( lBoldt& Ghosh 19991: iLevinson 2000l) . radio-loud 
AGN with jets seem privileged in that the most pronounced 
excess in arrival directions of UHECRs is positionally 
centered in the vicinity of the FR-I radio galaxy C en- 
taurus A ([Abraham et al. 20081; [Abraham etal. 2009h . '' 
FR-I radio galaxies, including their aUgned counterparts 
(BL Lac objects), radiate a volume- and time-averaged 
emissivity of » 10"^^- 10^^ erg Mpc~^ yr"' in nonther- 
mal 7 rays dPermer & Razzaque 20101) . and FR-I radio 
galaxies are found within the sa 100 Mpc Greisen- 
Zatsepin-Kuzmin radius. If comparable power goes 
into the acceleration of UHECRs, then BL Lac ob- 
jects and FR-I radio galaxies have more than sufficient 
emissivity to power the > 10 E eV UHECRs, which re- 
quire ^ 10"^ e rg Mpc"^ yr"' ( Waxman & B ahcall 1999t 
iBerezinskv et al. 2006; Murase & Takami 2009) . 

* See Takami & Sato 2009 for discussions on issues of UHECR anisotropy 
for protons, and Lemoine & Waxman 2009; Abreu et al. 2011 for heavy 
nuclei. Note also the potential contribution from the background Centaurus 
supercluster pointed out by Ghisellini et al. (2008). 



2.1. Two-Component Spectra of Blazars and Radio Galaxies 

In the framework of one-zone leptonic synchrotron/SSC 
models, the double-humped SED of blazars and radio galaxies 
(misaligned counterparts), plotted as log(z^/v) vs. log(i^), can 
be well described by a low-energy synchrotron component 
and a high-energy inverse-Compton (IC) curve. Each com- 
ponent is characterized by the peak synchrotron flux I'Fl at 
peak synchrotron photon energy e,, and the peak IC flux i/F^ 
at peak IC energy ec, respectively. The spectrum of radio- 
loud AGN is often highly variable, and rises and decays with 
variability time fvar- Here, we define the variability time as 
the shortest timescale in which the flux shows a significant 
factor-of-two change, and assume co-spatiality of the highly 
variable emissions in the different energy bands, as indicated 
by correlated variability dAckermann et al. 201 Ih . 

Here we consider such a one-zone model. The emis- 
sions are produced by electrons in a spherical blob moving 
relativistically in a jetted geometry: the high-energy emis- 
sion is produced via Compton scattering off the local syn- 
chrotron photons that are generated by the electrons in the 
blob. The synchrotron (IC) luminosities at peak energy are 
jj^O ^ A-Tid1{£s(C)Ff'"^) and the Compton dominance param- 
eter is defined as Aq = 1^11^% ~ {scFe^/iSsFP). The electron 
distribution is typically assumed to consist of power-law seg- 
ments. As long as both the vF^ and vF^ energy fluxes origi- 
nate from electrons with the same mean comoving-frame en- 
ergy 7^meC^, the magnetic energy density can be rewritten in 

terms of Ac as B'VStt - U^/{AttR''^5^cAc). Then the Doppler 
factor is given in this relativistic spherical blob formalism by 
the expression (ghisellini et al. 1996) 
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3(5(ec/£.s) 

provided the Compton scattering takes place in the Thomson 
regime, which applies when 5 >^ 5t = 2^/3\/es^c/n^c^{\ + 
z)^ Here the critical magnetic field is defined as Bq = 
mlc^/eh ~ 4.4 x lO'^ G (e.g., iBrainerd & Petrosian 1987|) . 
The above equations are derived by using the common 
relations £s/meC^ « S{B' /Bq)jI /(1 + z), and ec/nieC^ w 

(4/3)7^ £s (i.e., the typical fluid-frame Lorentz factor of 
electrons radiating near the peak synchrotron and SSC fre- 
quencies is 7! » ( ■\/3/2)-^/£c/£.; in the Thomson limit) (e.g., 
ISikora et al. 2009h . 

For high-peaked BL Lac objects, the Compton scattering 
often occurs in the KN regime, where more detailed modeling 
is required. Nevertheless, equations (1) and (2) demonstrate 
that source parameters such as & and B' can be determined 
from the double-humped SED.*^ 

^ The inequality is derived from 47^ej(l + z)/(5 < niiC- and ec ~ 

(4/3)7,;'£.,. 

** Other parameters such as the acceleration efficiency r} {t!^^^ = 
r]^'meC^ /(eB' c)) depend on details of the electron distribution. For exam- 
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Table 1 gives measured and inferred properties for blazars 
and radio galaxies with good multiwavelength coverage that 
are well described by a synchrotron/SSC model, where de- 
rived values of magnetic fields and Doppler and Lorentz fac- 
tors based on detailed synchrotron and SSC modeling are 
taken from the literature (rather than using equations 1 and 
2). In Table 2, we also show parameters derived with equa- 
tions (1) and (2) (only when S >6t i& satisfied). For blazars, 
we assume that F « (5, whereas values of the angle of the jet- 
ted emission with respect to the observer inferred from obser- 
vations are considered for radio galaxies. Also, R' « cf^^j. = 

C<^fvar/(1+Z)is USCd. 

The SED modeling has much improved thanks to the con- 
stantly increasing multi-wavelength coverage, which also al- 
lows simultaneous multi-band observations. Nevertheless, 
there is still some degree of scatter among parameter sets 
obtained by different groups, partially related to unavoid- 
able parameter-degeneracy with respect to the observables. 
In some cases, there is even significant scatter between de- 
rived parameters for the same state of a specific source (Cen 
A) or between different states of the same source (Mrk 501). 
In particular, the whole SED of Cen A up to the TeV-band 
cannot be fitted with a unique parameter set. Abdo et al. 
(2010d) show different models and parameter values for the 
same data, which reflect the limitations of fitting single epoch, 
single component SEDs for derivations of source parameters. 
However, our conclusions are not affected by uncertainties of 
the poorly constrained parameters, as we discuss in more de- 
tail below. 

2.2. Implications for UHECR Acceleration 

The Hillas condition dHillas 1984t) limits the maximum ac- 
celerated energy of ions with charge Z to 

E'Z-'^KZeB'TR', (3) 

in order that the particle Larmor radius is smaller than the 
characteristic size scale R' < ct'^^ = C(5fvar/(1 +z)- The in- 
equality is replaced by an equality in our estimates. Equa- 
tion (O can be rewritten using the Thomson-limit relations 
given above when (5 « F, in which case one gets 
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4(l;) 



.sn1/4. 1/2^1/2 



vai" Oq £s 



l,cV^Afef 



(4) 



Notice that Zs™^ is defined in the cosmic rest frame at UHECR 
production. 

In Table 2, the maximum proton energy is estimated using 
equation (3) and parameters given in Table 1, which in turn are 
based on the results of synchrotron/SSC model fits for these 
sources found in dedicated modeling papers (rather than using 
equations 1 and 2). Here, note that other losses and details of 
the acceleration process could limit the maximum particle en- 
ergy further For the cases considered in Table 1, it is barely 
possible to accelerate protons up to ^ 10 * eV, whereas Fe 
nuclei could easily each reach > 10 eV provided that they 
can survive photodisintegration. Alternately, > 10^*' eV pro- 
ton acceleration could occur transiently during rare bursts or 
flares CiPermer et al. 20091: iMurase & Takam i 2009), though 
according to Table 2 it might still be difficult even for bright 
flai-es from Mrk 501 and PKS 2155-304. As noticed above, 

pie, for the simple powe r-law injection of accele rated particles, large values 
of T] are often suggested flnoue & Takahara 1999) . which may not be the case 
for injection with multiple power-law segments. 



there is the large scatter in parameters due to uncertainties. 
However, even with the allowed spread in the parameter val- 
ues (see Table 2), this conclusion seems robust. 

A similar conclusion is also reached when considering lu- 
minosity requirement s for B L Lac objects and FR-I radio 
galaxies (i Dermer &Razzaque 2010t iGhisellini et al. 2010) . 
In Ghisellini et al. (2010), physical parameters were obtained 
via spectral modeling using their one-zone leptonic model of 
all blazars with known redshift detected by the Fermi satel- 
lite during its first 3 -month survey. The inferred magnetic 
luminosity of BL Lac objects is typically L^ ^ 10 b\ erg s"' 
and almost all of them satisfy Lb < 2 x 10 '' 8\ erg s"', where 
5\ = (5/10. On the other hand, the required magnetic lu- 
minosity for UHECR acceleration to lO-^L^o ^^ i^ Lb w 
2 X 10'^^F2(£™^)2z-2 erg s"' where Fi = F/10. Hence, it 

also suggests difficulties in acceleration of > 10 eV pro- 
tons in typical BL Lac objects, though the simple SSC model 
cannot be so simply applied to lower-peaked BL Lac objects 
where values of B' and 6 are not well-defined due possibly to 
external Compton scattering components. 

If BL Lac objects and FR -I radio galaxies, as has 
often been considered (e.g.. [Tin vakov & Tkachev 20011: 
iBerezinskv et al. 20021; iDermer eta l. 2009), accelerate the 
UHECRs, then an ultra-high-energy (UHE) proton origin 
of the highest-energy cosmic rays is disfavored from spec- 
tral modeling if the standard synchrotron/SSC model is cor- 
rect. Heavier nuclei can, however, be accelerated up to 
ultra-high energies. The composition of UHECRs is an 
open question, with both proton and heavy-ion dominated 
compositions having been claimed to be compatible with 
HiRes (Abbasi et al. 2010) and the Pierre Auger Observatory 
(PAO) ( Abraham et al. 201 0) data, respectively. As seen here, 
the standard model for 7-ray emission from BL Lac objects 
and FR-I radio galaxies suggests a transition from proton to 
heavy-ion dominated composition at^(10' -lO'') eV. 

We have assumed that 7 rays from TeV blazars and ra- 
dio galaxies are from leptonic Compton scattering in order 
that synchrotron theory can be used to derive the various pa- 
rameters. Thus the hadronic 7-ray flux must be consider- 
ably smaller for a consistent interpretation. Sufficiently high- 
energy protons and nuclei interact with synchrotron photons 
in the jet via the photomeson process, with photopion pro- 
duction efficiency fo-, for cosmic -ray protons estimated to 
be (e.g.. lMurase & Beacom 2010ab 



fp 



2.3 X 10" 



2.5 
1+a 



■^7,46'vai-,4''l 



IkeV 



where E'*^ 1.6 x 10 



16 



(5) 
eV (ejl keV) '(5?(l + z)"^ is the typ- 
ical energy of a proton that interacts with a photon with 
£5 (where the proton energy is here defined in the observer 
frame). Also, a is the photon index at energies below or above 
Ej. For a w 1 .5, which is typical of BL Lac objects at Ep > E^, 

the photomeson production efficiency at ^ 10 eV becomes 
of order /p^ ^ 6 x 10" , which suggests that the photomeson 
process is inefficient for this kind of blazar (though it could 
be more efficient for low-peaked BL Lac objects and FSRQs). 
The efficiency can also be higher if F is lower, provided that F 
is consistent with synchrotron/SSC model fits and minimum 
Lorentz factor estimates inferred from, e.g., 77 opacity argu- 
ments. 
Roughly, half of the pions produced by photome- 



TABLE 1 

Measured and inferred properties of VHE blazars and radio galaxies 



ID 


Source 


z 


Epoch 


fvar 


gia) 


obs 


<'•' 


e.S"' 


!^.,f;''''[10"'°] 


R'<"'[10'=] 


B'W 


4°' 


^^^^[lo-iO] 


Ref. 










[s] 




/[deg] 




[meC^] 


[ergcm"^ s"'] 


[cm] 


[G] 


[nicC ] 


[erg cm"' s"'] 




I 


CenA(core) 


0.00183 


2009 


< 1.0 X 10^ 


1.0-3.9 


2.0-7.0/15-30 


(0.8-400) X 10-' 


(0.8-4000) X 10"' 


0.09-4.5 


3.0-11.0 


0.02-6.2 


0.17-(8.3x 10^) 


0.025-8.5 


1 


2 


M87 


0.00436 


2009 


1.7 X 10^ 


3.9 


2.3/10 


4x 10^ 


1.6 X 10"' 


0.06 


14.0 


0.055 


18.6 


0.068 


2 


3 


NGCI275 


0.0179 


Oct. 2010'' 


8.6 X 10" 


2.3 


1.8/25 


960 


2.4 X 10"^ 


0.9 


2x 10^ 


0.05 


2.9 X 10^ 


0.3 


3 


4 


NGC6251 


0.024 


- 


- 


2.4 


2.4/25 


2x lO" 


6.5 X 10"' 


0.012 


120 


0.037 


7.3 


0.047 


4 


5 


Mrk421 


0.03 


19 March 2001 


1.0x10^ 


80 


80 


9.3x10" 


0.005 


7.4 


3.0 


0.048 


8.1 X 10" 


7.0 


5 


6 


Mrk501(h.<'',1997) 


0.0337 


16 April 1997 


7x 10^ 


14-20 


14-20 


(7-300) X 10" 


0.3-0.5 


8.0-8.5 


1.0-5.0 


0.15-0.8 


(1.4-2.6) X 10" 


2.9-3.4 


6,7,8 


7 


MrkSOl (!.''', 1997) 


0.0337 


7 April 1997 


- 


15 


15 


6x 10^ 


0.002 


0.63 


5.0 


0.8 


4.4 X 10' 


0.4 


6 


8 


MrkSOl (l.''',2007) 


0.0337 


2007 


- 


25 


25 


1 X 10^ 


0.002 


0.63 


1.0 


0.31 


4.4 X 10' 


0.4 


9 


9 


MrkSOl (!.''' ,2009) 


0.0337 


2009 


3.5 X 10^ 


12-25 


12-25 


(6-90) X lO" 


0.002 


0.55-0.63 


1.0-130 


0.015-0.34 


(1.3-4.4) X lO' 


0.3-0.4 


7,10,11 


10 


lESI959+650(li.'"') 


0.047 


Sept2001-May2002 


(2.2-7.2) X lO" 


18-20 


14-20 


4-5 X lO" 


(0.07-8) X 10"^ 


1.0-3 


5.8-9 


0.04-0.9 


8 X 10'"" 


0.2-2 


12,13 


11 


1ESI959+650(1.<'') 


0.047 


23-25 May2006 


8.64 X lO" 


18 


18 


5.7 X lO" 


0.003 


2.6 


7.3 


0.25-0.4 


1.2 X lO' 


0.22 


14,15 


12 


PKS2200+420/BL Lac 


0.069 


- 


- 


15 


15 


900.0 


5.3 X 10"' 


0.76 


2.0 


1.4 


1.6 


0.4 


14 


13 


PKS2005-489 


0.071 


- 


- 


22 


22 


1.3 X 10" 


4.7 X 10"^ 


1.5 


8.0 


0.7 


3.6 X 10^ 


0.07 


14 


14 


WComae 


0.102 


7-8 June 2008 


5400 


20-25 


20-25 


(1.5-20) X lO" 


8.0 X 10"^ 


0.4 


3.0 


0.24-0.3 


8.1 X 10^ 


0.15 


14,16 


15 


PKS2I55-304 


0.116 


28-30 July 2006 


300 


110 


110 


4.3 X lO" 


4x10"" 


2.13 


0.86 


0.1 


9.7 X lO' 


20.0 


5 



'"': parameter value from the SED modeling in literature (see references); *''': for blazar sources <5 Ri F and 9j Ri l/F; *"': high (h.) and low (1.) state; References: 1- Abdo et al.(2010d) (see Figure 5 and Table 2 in 
the paper for the different models), 2- Abdo et al. (2009c), 3- Abdo et al. (2009b) (see also Brown & Adams 2011), 4- Migliori et al. (2011), 5- Finke et al. (2008), 6- Plan et al. (1998), 7- Acciari et al. (2011), 8- 
Katarzynski et al. (2001), 9- Albert et al. (2007b), 10- Anderhub et al. (2009), 11- Abdo et al. (2011c), 12- Tagliaferri et al. (2003), 13- Krawczynski et al. (2004), 14- Tavecchio et al. (2010), 15- Tagliaferri et al. 
(2008), 16- Acciari et al. (2009b). 



CI 



BL Lac Objects as UHECR Sources 



TABLE 2 
Inferred properties of VHE blazars and radio galaxies 



ID 


Source 


di 


il,[10*] 


L^[IO"] 


Ac 


II, 


sp 


S"> 


^'""[10"] 


B'"" 


£r"'/Z[10"'] 






[Mpc] 


[erg cm"- s"'] 


[erg cm"- s"'] 










[cm] 


[G] 


[eV] 


1 


CenA(coi-e) 


3.7 


(0.15-7.3) X 10"^ 


(0.04- 14) X 10"* 


0.3-1.9 


890-2.1x10* 


9.9-(6 X 10"*) 


0.12-3.7 


3.0-12 


0.02-9.1 


0.01-4 


2 


M87 


16.7 


2.0 X 10"* 


2.3 X 10"* 


1.1 


9.3 x 10' 


0.006 


2.7 


20 


0.021 


0.05 


3 


NGC1275 


75.3 


0.06 


0.02 


0.35 


960 


0.005 


- 


- 


- 


5 


4 


NGC625I 


104 


2 X ur^ 


6.6 X 10"-' 


3.3 


2.9 x 10' 


0.007 


- 


- 


- 


0.3 


5 


Mi-k42I 


130.0 


1.5 


1.4 


0.95 


3.4 X 10-' 


74 


- 


- 


- 


0.3 


6 


MrkSOI (h.<«',1997) 


146.0 


2.0-2.2 


0.7-0.9 


0.36-0.41 


(1.4-2.5) X 10' 


(3.0-3.5) X 10' 


- 


- 


- 


0.1-2 


7 


Mrk50I(l>\l997) 


146.0 


0.2-0.4 


0.1-0.2 


0.44-0.63 


(0.08-I.3)xI0* 


100-1700 


- 


- 


- 


2 


8 


MrkSOl (1>',2007) 


146.0 


0.2 


O.I 


0.63 


1.3x10* 


100 


- 


- 


- 


0.2 


9 


Mi-kSOI (l.<«',2009) 


146.0 


0.1-0.2 


0.08-0.1 


0.55-0.63 


(0.7- 1.3) X 10* 


58-100 


- 


- 


- 


0.2-0.7 


10 


lES1959+650(h.'«') 


206 


0.5-1.5 


O.I-I.I 


0.2-0.8 


(2.7-9.5) X 10* 


27-910 


- 


- 


- 


0.1-3 


11 


lES1959+650(l.'-5') 


206 


1.3 


O.I 


0.08 


6600 


66 


- 


- 


- 


1-2 


12 


PKS2200+420/BL Lac 


307.0 


0.8 


0.45 


0.53 


2.8x10' 


0.006 


- 


- 


- 


I 


13 


PKS2005-489 


316.0 


1.8 


0.07 


0.04 


7.6x10' 


1.5 


- 


- 


- 


4 


14 


WComae 


464.0 


1.0 


0.38 


0.38 


8.7x10' 


3.1 


7.2 


3.0-3.7 


2.1-2.6 


0.4-0.7 


15 


PKS2 1 55-304 


533.0 


7.2 


68 


9.4 


1.3x10* 


24 


- 


- 


- 


0.3 



'"*: obtained from equation (2); '*': St = 2\/3\/£c£s/'"?c''(l +z)', *^': obtained from equation (1); '''': calculated assuming R' : 
7^ Si si I2yjeclss', "^': obtained from equation (3) using B' , F, and R' reported in Table 1; '*': high (h.) and low (1.) state. 



;cr^„ = c5rvar/(H-z); 



equivalent) neutrino luminosity at given Ey f 



is estimated to be EU^ 
lO'*^ erg s" 



son production are charged, and each neutrino car- 
ries ^ 1/4 of the pion energy, so the total (isotropic- 

0.05£,, 
- 3.8 X 
' fp^-2(El^ jVd ^ erg s"') for a source satisfy- 
ing fp^ < 1 like BL Lac objects, where fp^ = 10" fp^-i is 
used. Then, one finds that the neutrino flux from an indi- 
vidual source is typically too low to be detected with Ice- 
Cube. One can also see that the cumulative background flux 
from high-peaked BL Lac objects is low. The UHECR en- 
ergy input in the local universe is ~ 5 x 10 erg Mpc~^ yr"' 
at lO'^ eV (iMurase & Takami 2009h . so that assuming that 
such BL Lac objects and FR-I galaxies are the main UHECR 
sources, the expected cumulative muon neutrino background 
flux is estimated to be (iMurase & Beacom 2010al) 



£>.^ 



10"'° GeV cm"^ s" 



fjn(20E,) 
10"^ 



p2-p r 



(6) 

where p is the cosmic -ray spectral index and /, is a pre- 
factor coming from the redshift evolution of the sources. Low 
photomeson production efficiencies also follow if the UHE- 
CRs are heavy nuclei, whose losses are dominated by pho- 
todisintegration (see below). More luminous blazars, includ- 
ing low-peaked BL Lac objects, may however lead to higher 
photomeson production efficiencies, so that the cumulative 
neutrino backg round cou ld be dominated by this class of 
AGN ( Mucke et al. 2003h . 

Our evaluation is based on the standard synchrotron/SSC 
model for BL Lac objects and FR-I radio galaxies. One 
could abandon the standard synchrotron/SSC model and con- 
sider a highly magnetized, ^ 10-100 G jet model, which 
is needed in hadronic blazar models to accelerate protons 
to > 10^° eV (lAharonian 20001) . Correspondingly, the min- 
imum magnetic luminosity is estimated to be L^ w 2 x 
' r?(£'l™A) , which is larger than the typical syn- 



lO''^ erg s 



^z/ rma 
l'-^p,2 



lO"*^ erg s" 



In 



chrotron luminosity of BL Lac objects, U 
the hadronic models, 7-ray emission is attributed to proton 
synchrotron radiation and/or proton-induced cascade emis- 
sion, which leads to the requirement that the UHECR lumi- 



nosity is Z^uHEC R <', L^ = ^cLj. In the proton synchrotron 
blazar model (Aharonian 2000"; ' Miicke & P rotheroe 20011: 
Mucke et al. 2003), which is typically viable for high-peaked 
BL Lac objects, proton synchrotron radiation is emitted up to 
energies e^ « 4(5i/(l +z) TeV (in the limit that the maximum 
energy is determined by the synchrotron cooling) in efficient 
Fermi acceleration scenarios. The photomeson production ef- 
ficiency for protons is strongly dependent on S, but the con- 
dition Lb ^ EL suggests that the synchrotron energy-loss is 
dominant at ultra-high energies where the proton synchrotron 
radiation is typically prominent at ^ TeV energies (though 
the photohadronic cascade component may become relevant 
at lower energies). The strong magnetic field also suppresses 
electronic SSC emission because fewer electrons are needed 
to generate the same synchrotron flux. 

2.3. Survival of Nuclei in the Source 

If the standard synchrotron/SSC scenario holds for TeV 
blazars and their misaligned counterparts, then the protons 
can hardly reach 10^° eV, as shown in Table 2 (see E™"" ob- 
tained by detailed modeling in the literature). For BL Lac 
objects and FR-I galaxies to be the steady sources of UHE- 
CRs, therefore, UHECRs would primarily be heavier nuclei. 
In such a scenario, one has to exami ne whether ions can 
survive photodisintegration losses (cf. IMurase et al. 2008al : 
'Wang et al. 20081: iPe'er et al. 20091) . The photodisintegration 
opacity is estimated similarly to the photomeson production 
efficiency. Approximating the photodisintegration cross sec- 
tion by the giant dipole resonance (GDR) cross section as 
<ta-i ^ crGDR^(e - £GDR)AgGDR , for a sufficiently soft pho ton 
spectrum (with a > 1), we get (IMurase & Beacom 2010aL see 
also Murase et al. 2008a for the non-GDR effect) 



7A7 ■ 



'A7 



2crGDR AegDR 



1-1- a EgDR 47r^'*fvarC^£* V^^' 



(7) 



where f L is the photodisintegration interaction time, ctgdr ~ 
1.45 x 10"^^A cm^, egtik ~ 42.65A-°2i MeV (for A > 4), 
AecDR - 8 MeV, and E^ w 0. 5 S\m ac^ s cDR/eK I +zT^ (in 



Murase et al. 



the observer frame). Then we numerically find 

, 2.5 \ . . / p- N-i 



l+a/'^''^'"^'™" 



''■4'^iATk^) 



a-l 



(8) 
where E'^ ~ 4.8 x lO'* eV (A/56f™(£,/l keV) ^5](\+zr^ 
is the energy of a nucleus that typically interacts with a pho- 
ton with £j. Hence, heavy nuclei with Ea ^ (Z/26)10 eV 
(given in the observer frame) undergo some photodis- 
integration reactions unless 6 is high enough. The 
nucleus survival condition TA-^iEA) < 1 gives 6 > 

ll{Z/26f\A/56)-'-''\u^,,)'lXl!l{eJl keV)-«'(l ^zT'" 
(for a^ \ .5), but significant photodisintegration loss is easily 
avoided for reasonably large bulk outflow Doppler factors. 

Recalling from equation (|5]l that the photomeson produc- 
tion efficiency has the same dependence on 6, we can con- 
clude that when heavy nuclei survive photodisintegration, 
the photomeson production efficiency is so low that the 
correspond ing neutrino and 7-ra y fluxes are not easily de- 
tected (Mu rase & Beacom 2010al but see also Murase & Bea- 
com 2010b). 

3. EXTREME TEV BLAZARS AND INTERGALACTIC CASCADES 

The 77 opacity argument allows us to place constraints on 
5 in BL Lac objects obser ved at TeV e nergies, and requires 
(5 > 60 for PKS 2155-304 (iBegelman et al. 2008) for the ma- 
jor July /August 2006 TeV flares dAharonian et al. 2007ah . and 
5 > 100 to be furthermore consistent with synchrotron/SSC 
model fitting for di fferent models o f the extragalactic back- 
ground light (EBL; F inke et al. 2008 ). Another important fact 
is that VHE photons can interact with the cosmic photon back- 
grounds. VHE 7 rays produce electron-positron pairs via 
77 pair creation, and the resulting high-energy pairs make 
high-energy photons via Compton scattering. Hence, the 
cascaded 7 rays, which are often called pair echoes (e.g., 
Plaga 1995; Murase e t al. 2008b) and/or pair ha loes (e.g.. 



Aharonian et al. 1994t iNeronov & Semikoz 2007h . are ex- 
pected at GeV-TeV energies. In particular, 7 rays with en- 
ergies below ^100 TeV are likely to leave structured re- 
gions of the universe, inducing the cascade in the void re- 
gion (Murase et al. 2008b). 

Many 7-ray blazars show variability, and often dis- 
play spectacular flares. Some of them are ultra-variable, 
as se en in multi-TeV flari n g episodes from PKS 2155- 
304 ([Aharonian et al. 20074 lAhar onian et al. 2009b') , Mrk 
501 (E-lbert et al. 2007b) and Mrk 421 (Galante et al. 201 lb . 
Such rapidly varying 7-ray emission should be produced in 
the blazar region. This is because the IGMF will introduce a 
significant time spread in the cascade radiation, which seems 
incompatible with rapidly varying emissions. For example, 
consider primary 10 TeV 7 rays, so the Compton-upscattered 
cosmic microwave background (CMB) photons have E^ w 



(4/3)7'2ecMB 
,1/2 



8 GeV j'j. Based on the lower limits 

of B^cvXtLl > 10" ' ^- 10"'^ G Mpc'/^ (iDolagetal. 201 U 
iDermeretal. 2OIII: iTakahashi et al. 2012h obtained for IBS 
0229H-200 (where the term Acoh is the coherence length of 
the magnetic field, and this relation is understood to apply 
when Acoh is smaller than the cooling length for GeV produc- 
tio n), the timescale of the ^0.1 TeV pair ech o is estimated to 
be (iMurase et al. 2008bl; iDermer et al. 201 lb . 

AfiGV-^ ...-^-2 02 .X /,nn^^„„^., , .^-I 

where A 

the void IGMF Bigv is defined in the frame of the Hubble 



1 .4 yr £-2 1^2^^ _,7(A^V100 Mpc)(l +zT' , (9) 
^j is the mean free path for the 77 pair creation, and 



flow. Therefore, the highly variable VHE radiation from these 
BL Lac objects is likely to be either leptonic synchrotron or 
SSC, or proton synchrotron radiation produced in the jet (even 
though another slowly variable component may be produced 
by the secondary emission). 

This conclusion does not however hold for a fraction of 
blazars and radio galaxies from which prominent variabil- 
ity has not been seen. An interesting source is the ex- 
ti-eme TeV blazar lES 0229H-200, which has a hard VHE 
component extending to > 10 TeV, but has not been re- 
ported to b e variable in observations taken over a pe riod 
of 3 - 4 yr (lAharonian et al. 2007bt iPerkins et al. 2"oTol) . If 
the apparent absence of the variability comes from observa- 
tional limitations and if fast variability is present, the emis- 
sion should be produced in/near the blazar region. If it is 
the case that there is no rapid variabilit y, the observed com- 
ponent may come from an extended jet CBottcher et al. 20081 
see also Section 4 for further discussion on the 7-ray emis- 
sion region). In addition, as we see in this paper, a 
slowly variable component can be 7-ray-induced intergalac- 
tic cascade emission. Furthermore, if it is non-variable, 
proton-induced intergalactic cascade emission (i.e., inter- 
galactic cascades caused by UHE 7 rays and pairs gen- 
erated via the photomeson production with the CMB and 
EBL) can be responsible for the observed emission. These 
cascade emissions may confuse interpretation of the mini- 
mum bulk Lorent z factor from 77 opacity arguments and the 
level of the EBL (lEssev & Kusenko 20101; lEssev et al. 20101: 
lEssey etal. 201 Ih . The intergalactic cascade components 
could be present not only in extreme blazars, but also on 
longer timescales in other blazars and radio galaxies. For ex- 
ample, a slowly variable (~ month) emiss ion at GeV energies 
was observed from Mrk 501 (lAbdo et al. 201 Ic) . which could 
arise from the i ntergalactic cascade i nduced by variable TeV 
source photons (i Neronov et al. 2011]) . 

Here we focus on the intergalactic cascade scenarios in 
order to explain hard VHE spectra of extreme TeV blazars 
whose variability is apparently absent. We calculate the cas- 
cade emission by solving the Boltzmann equations, where 77 
pair creation, IC scattering, synchrotron radiation, a nd adia- 
batic energy loss are taken into account (iLee 1998h . As for 
proton propagation, we directly solve the equation of motion 
of protons one by one, with photom eson production simu- 
lated by SOPHIA (iMiicke et al. 20001) and the Bethe-Heitler 
process in cluded to treat interactio ns with the ambient pho- 
ton field dChodorowski et al. 1992h . Then the electromag- 
netic cascade is calculated separately. Fo r the EBL model, we 
employ the low-IR an d best-fit models (iKneiske et al. 20()4i 
iKneiske & Dole 20101 see Finke et al. 201 1 for detailed dis- 
cussions on the EBL). We focus on the possibility that the 
cascade interpretation is a viable explanation of VHE 7-ray 
spectra of extreme TeV blazars. Here, the IGMF in voids 

has to be weak enough (BigwXIH < 10"'^ G Mpc'/-) that 
the cascade radiation at TeV energies is not suppressed by 
the IGMF if it is to make the measured flux in the VHE 
range. On the other hand, the void IGMF cannot be be- 
low ^ 10" G Mpc''2 due to constraints from Fermi (e.g. , 
'Dolag et al. 201 1*; 'Permer et al. 20 llt iTavecchio et al. 201 ll: 
Taylor et al. 201 1; Ahlers & SalvadoJOni). 

An important point is that cosmic magnetic fields are al- 
most certainly inhomogeneous. Whereas one may expect very 

weak IGMFs in voids, BicyXl'^ < 10"'' G Mpc'''^ the struc- 
tured region of the universe is likely to be significantly mag- 



BL Lac Objects as UHECR Sources 



netized. Cl usters of gala xies are known to have Beg ^ 0.1 - 
1 /LtG (e.g., IVallee 20041) . and recent simul ations have sug- 
gested that filam ents have Beg ~ 1 - 10 nG dRvu et al. 20081: 
iDas et al. 20081 and see also, e.g., Donnert et al. 2009), which 
are larger than levels expected for the IGMF in voids, and 
galaxies including AGN are likely located in these structured 
regions of the universe. The mean free path of < 100 TeV 
and > 3 EeV 7 rays is larger than ^ Mpc (e.g.. lDermer 20071) . 
so that one may expect that the cascade emission induced by 
VHE/UHE primary 7 rays is primarily developed in the voids. 
On the other hand, ions must propagate in the clusters and/or 
filaments, so that they are deflected (and delayed) by their 
magnetic fields. Indeed, as demonstrated by a number of au- 
thors, the structured EG MFs play a crucial role on propaga - 
tion of UHECRs (e.g., ITakami et al. 20061: iDas et al. 2008h . 
and this is even more so the case for lower-energy cosmic 
rays. 

3.1. Cascades by Primary VHE/UHE Gamma Rays 

SEDs of high-peaked BL Lac objects are generally 
well reproduced by the standard one-zone electronic syn- 
chrotron/SSC model. Among them, extreme TeV blazars 
have the hardest VHE 7-ray spectra at ^ 1-10 TeV ener- 
gies, as indicated by deabsorption of the measured 7-ray spec- 
trum based on conventional EBL models (discussed below) 
and supported by non-detections of GeV 7 rays by Fermi. 
Also, in some cases (RGB J0152H-017, lES 0229H-200 and 
lES 0548-322), the optical/UV data show a rather steep spec- 
tru m which is thought to b e the emission from the host galax- 
ies (iTavecchio et al. 20111) . Although the synchrotron compo- 
nent of extreme TeV blazars seem unremarkable at the opti- 
cal/UV band, in these cases, comparison between optical/UV 
and X-ray data requires a strong roll-off of the nonthermal 
spectrum below the X-ray band, sugge sting that F^ozv^l^ for 
lES 0229H-200 (ITavecchio et al. 201 Ih . 

It is possible to explain such hard 7-ray spectra by the SSC 
model, but extreme parameters seem necessary compared to 
cases of typical, variable high-peaked BL Lac objects. It of- 
ten suggests a very narrow-range energy distribution of elec- 
trons, and unusually large values of (5 ^ 10 - 10^ may be 
necessary to avoid the Klein-Nishina suppression. For lES 
0229H-200, PKS 0548-322 and lES 0347-121, extreme values 
of the electron minimum Lorentz factor of 7(. „, ^ 10-10^ are 
required from spectral modeling, where the hard SSC spec- 
trum, Fe oc E'"^! ^ (in the Thomson regi me) can be expected in 
the VHE range (ITavecchio et al. 201 Ih . 

There are several alternate blazar models that predict very 
hard ~ 30-100 TeV 7-ray emission. In the hadronic 
model, the proton synchrotron process leads to multi- 
TeV emission if the outflow is ultra-relativistic, F ~' 10 - 
10 . Then, further hardening may be caused by inter- 
nal a bsorption due to some sof t photon field outside the 
blob (IZacharopoulou et al. 201 Ih . Another possibility to 
make hard TeV emission is electromagnetic radiation pro- 
duced by nonthermal elec trons in the vacu um gap of the black 
hole magnetosphere (e.g., Levinso n 20001) . 

Bottcher et al. (2008) suggested that hard VHE emission 
originates from CMB photons Compton-upscattered by rel- 
ativistic electrons that are accelerated in the extended jet. 
In this model, if the electron spectrum is hard, p ^ 1.5, 
the resulting number spectrum of the IC emission has Fe oc 
£'-(i+p)/2 ^ £-5/4^ which is compatible with the observed 
VHE 7-ray spectrum. The same process could be important 
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Fig. 1. — Spectra of VHE 7-ray-induced cascade emission for various 
source redshifts. We assume the total 7-ray luminosity of L-y = Iff*^ erg s~' 
with P = 2/3 and £""« = 100 TeV. The low-IR EBL model of Kneiske et al. 
(2004) is used here. 
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Fig. 2. — Spectra of VHE 7-ray-induced cascade emission for various 
intrinsic photon spectra. The source redshift is set to z = 0.14. 



for recoUimation sho cks (e.g., iBromberg & Levinson 20091 : 
iNalewaiko & Sikora2 009), or acceleration at knots and 
hotspots, noting that variations on much longer timescales can 
be expected in these models. 

When VHE 7 rays are emitted from a source, they in- 
duce an electromagnetic cascade in intergalactic space. This 
cascade unavoidably accompanies spectral production of ex- 
treme TeV blazars as long as the IGMF in voids is weak 
enough. To demonstrate this, we show in Figures [T] and |2] 
the VHE 7-ray-induced cascade emission for sources at var- 
ious redshifts. In Figure [H primary source photons with 
Fe oc E-P with /3 = 2/3 and £"" = 100 TeV (in the source 
rest frame) are assumed. One sees that the observed cutoff 
due to the EBL becomes lower for more distant sources since 
the 77 pair-creation opacity increases. In Figure|2] a different 
photon index ifi = 5/4) and/or a different maximum energy 
(-^■max _ jqI 5 ^gy-j ^j-g assumed for comparison, which causes 
slight differences in spectra. As indicated by equation (9), the 
intergalactic cascade emission induced by primary 7 rays will 
be slowly variable or almost steady. 

Another way to have 7-ray induced emission involves UHE 
7 rays produced in blazar jets or radio galaxies. Such a 
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Fig. 3. — Spectra of UHE 7-ray-induced cascade emission for various 
source redshifts. We assume L^y = 10^^ erg s"' at lO'**^^ — lO"^^ eV. 



case is shown in Figure 3 assuming much higher injected 
photon energies than before. Here we assume an injection 
spectrum centered at 10 EeV spanning one decade. The 
photomeson production by UHE protons, which can be ex- 
pected in hadronic models, leads to UHE photons with energy 
E^ « O.lEp ~ 10'' eV Ep^20- In the synchrotron source in 
which UHE protons are accelerated, one may expect that the 
synchrotron self-absorption cutoff curtails the number of low 
energy phot ons impeding UHE photon escape from the emis- 
sion region (Murase 2009). A caveat of this model in our case 
is that generation of UHE 7 rays in the source requires accel- 
eration of UHE protons and moderately efficient photomeson 
production. As noted before, the photomeson production in 
the source may not be too efficient in high-peaked BL Lac 
objects, which implies that the required UHECR luminosity 
has to be very large. As can be seen, there is some notable 
differences at low redshifts z < 0.1 due to the longer effective 
energy loss length of UHE 7 rays, but the received spectra are 
not strongly sensitive to the energy at which the photons are 
injected for higher redshift sources. 

Note that the intergalactic cascade scenario makes a non- 
variable or slowly variable component, even when the 7- 
ray emission made in the jet contributes to a separate 
highly variable component. Although there is no strong ev- 
idence of time variability for several extreme TeV blazars, 
future se nsitive observat ions by CTA (A ctis et al. 2 011), 
HAWC iSandoval et al. 2 009). LHAASO (ICao 2010 ). or 
SCORE (iHam pf et al. 201 1.) will be crucial for identifying a 
slowly variable 7-ray emission component. 

3.2. Cascades by Primary UHECRs 

In the previous subsection, we considered cascade emission 
induced by primary 7 rays. VHE 7 rays at < 100 TeV, or UHE 
7 rays with energies > 3 EeV, where the opacity of the back- 
ground radiation is not so large, can leave structured regions 
of the universe, whereas cosmic rays should feel structured 
EGMFs in clusters and filaments. The deflection of cosmic 
rays by the structured region with size / ^ Mpc, magnetic 
field of Beg ^ 10 nG, and coherence length of Acoh ^0.1 Mpc 
(which may be typical of filaments; Ryu et al. 2008) is esti- 
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Therefore, the deflection by the structured EGMFs is not neg- 
ligible for cosmic rays with energies < 10 eV, since the de- 
flection angle is larger than the typical jet opening angle of 
9j ^0.1 ~' 6°. The corresponding time spread due to a struc- 
tured EGMF around the source (that is comparable to the time 
delay) is expected to be 
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which is unavoidable as long as cosmic rays pass through the 
structured region around the source, and it implies that the 
resulting cascade emission is essentially regarded as steady 
emission. Notice that the total time spread ATqr could gen- 
erally be longer than AfcR due to additional time spread by 
intervening structured EGMFs and the void IGMF 

In order to model the structured EGMFs, we have as- 
sumed a simplified two -zone model with struct ured EGMF 
and IGMF in voids (see'Taka mi & Murase 201 U for details). 
We model a cluster of galaxies by a sphere with the radius 

of 3 Mpc, and BEcir) = B()(l -t-r/r^)""'^, with Bo=l ^iG and 
re = 378 kpc. The magnetic field direction is assumed to be 
turbulent with the Kolmogorov spectrum and the maximum 
length of Amax = 100 kpc. In addition to the EBL, the in- 
frared background in the cluster is considered as the super- 
position of the S EDs of 100 giant elliptical galaxies calcu- 
lated by GRASIL (ISilva et al. 1998^, using fitting formula for 
the gas distribution (Rordorf et al. 2004). Fil aments are mod- 
eled by a cylinder with a radius of 2 Mpc ( Ryu et al. 20081) 
and a height of 25 Mpc. The magnetic field is assumed to 
be turbulent, which is described by the Kolmogorov spectrum 
with Beg =10 nG and A^ax = 100 kpc, although these val- 
ues are very uncertain. Some numerical simulations imply a 
large-scale coherent component of the magnetic field in fila- 
ments (e.g., Briiggen et al. 2005), which may deflect cosmic- 
ray trajectories even more effectively. UHECRs are injected 
from the center of the filament toward a direction perpendic- 
ular to the cylindrical axis in order to examine a relatively 
conservative case. Throughout this work, the IGMF in voids 
is assumed to be weak enough to be less important for cosmic- 
ray deflections. 

In Figure |4] we show our numerical results for the case 
£max _ 2q19 ^y gxpected in the standard synchrotron/SSC 
model of typical, variable BL Lac objects and FR-I galax- 
ies. One sees that the structured EGMFs play an important 
role by suppressing the resulting 7-ray flux by more than one 
order of magnitude compared to the case without them. In 
Figure |5] we show the case of Zs™" = 10 ' eV, which can 
be achieved in the hadronic model. While the Bethe-Heitler 
pair-creation process dominantly provides an electromagnetic 
component in Figure IH contribution of photomeson produc- 
tion is more important in Figure 5. In the filament case, the 
deflection angle of UHECRs around 10 eV is still less than 
the jet opening angle, so that the 7-ray flux is diluted by only 
a small factor. On the other hand, in the cluster case, be- 
cause UHECRs cannot be beamed, the 7-ray flux becomes 
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Fig. 4. — Effects of the structured EGMF on the 7-ray flux. We assume 
iuHECR = lO"*' erg s-', with E™" = lO" eV and p = 2. Here, as in the 
results on cascade emission induced by primary 7 rays, we use the isotropic- 
equivalent cosmic-ray luminosity at the source (defined for UHECRs above 
10 eV), which is related to the absolute (beaming-corrected) cosmic-ray 
luminosity, Luhecr,;, as Luhecr = (1 -cosSjO^'Luhecr.;- Here the assumed 
jet opening angle is 6^ = 0.1. The source redshift is set to z = 0.5. 
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Fig. 5.— Same as Figure 4, but with E™" = 10"" eV. 
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Fig. 6. — Effects of the structured EGMFs on the deflection of UHE 
protons. Relative contributions represent how much the apparent cosmic- 
luminosity at which cosmic rays enter the void region is diluted from El^ 
at the source. Note that a two-sided jet is considered throughout this work. 
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Fig. 7. — Spectra of UHE proton-induced cascade emission for various 
source redshifts. We assume Luhecr = 10 erg s"' with £"""' = 10 eV and 
p = 2. The source is assumed to be located in the filament with Beg = 10 nO 
and Aniax =0.1 Mpc. The low-IR EBL model is here assumed. 



almost isotropic and the corresponding flux is reduced ac- 
cording to the jet beaming factor (1 -cos9j) ~ 1/200 for 
9j = 0.1. The effects of the structured EGMFs are illustrated 
in Figure 6, where the relative contributions are calculated 
from two-dimensional Gaussian fits. Note that if we express 
the isotropic-equivalent cosmic-ray luminosity where cosmic 
rays leave the structured region as EL'^, then the relative 

contributions are {\-cosOj){EL^)/(EL^ '■'). In the filament 
case, isotropization becomes significant at ^ 10 eV rather 
than at ~' 10 ' eV for the cluster case. 

In Figure |7] we show resulting 7-ray spectra for various 
redshifts. Owing to the Bethe-Heitler process with energy- 
loss length ^ Gpc, UHE protons continue to supply electron- 
positron pairs for a longer distance than the photomeson en- 
ergy loss length of ^ 100 Mpc. As a result, the dependence 
of the proton-induced 7-ray fluxes on distance is much gen- 
tler than 7-ray-induced fluxes. Indeed, one sees that the rel- 
ative importance of the proton-induced 7-ray flux to the 7- 
ray induced flux increases with distance (compare Figure 7 



with Figures 1 and 3). Importantly for distant sources, the 
proton-induced cascade spectrum is much harder than the 7- 
ray induced spectrum, especially above TeV energies. Future 
VHE observations by CTA and HAWC are important to iden- 
tify the origin of UHECRs through detection of high-energy 
7 rays, as we now demonstrate for lES 0229-1-200 in the next 
subsection. 

In this work, we are interested in cases where IC cascade 
emission in voids is important in the VHE range, since it 
can explain hard VHE spectra of extreme TeV blazars as sug- 
gested by Essey et al. (2010). When pairs are mainly supplied 
via the Bethe-Heitler process, the timescale of secondary pho- 
tons produced by a proton beam roughly becomes 



At 



IGV ■ 



14yr£-2iB 



7,I1"IGV,-I7 



(Abh/Gpc)(1+z)- 



(12) 



which is more relevant than ATcr when the void IGMF 
is so strong that ATcr < Aficv is satisfied. Here, 
Abh is the Bethe-Heitler energy loss length. One 
should also keep in mind that the proton-induced GeV- 
TeV synchrotron emission from the structured region it- 
self, where the EGMFs are stronger, should also be ex- 
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pected (see iGabici & Aharonian 2001 iKotera et al. 20091: 
iKotera et al. 201 U and references therein). For a weak IGMF 
that is of interest in this work, its relative importance is some- 
what smaller when the volume filling fraction of the magne- 
tized region is taken into account. 

We have demonstrated the likely importance of the struc- 
tured EGMFs for proton-induced intergalactic cascade emis- 
sion. They are also important for UHE nuclei. Since nuclei 
with energy ZEp have the same deflection angle as protons 
with energy Ep, our results indicate that Fe nuclei should be 
significantly isotropized for all observed UHECR energies. 
For UHE nuclei, the photodisintegration energy loss length 
is ^ 100 Mpc, for which the energy fraction carried by 7 
rays and neutrinos is small as long as Zs™" is not too high. 
On the other hand, UHE nuclei supply high-energy pairs via 
the Bethe-Heitler process, whose effective cross section is 
'«bh,aO'bh,a ^ kbh,/;0'bh,/)(Z^M), which induces cascades in 
the same manner as UHE protons. Therefore, the intergalac- 
tic cascade signal, which is generated outside the source, is 
also important for sources of primary UHE nuclei. ^ 

3.3. Implications for TeV-PeV Observations 

In a wide range of EBL models, deabsorption of measured 
TeV blazar spectra leads to hard excesses at >TeV energies in, 
e.g., lES 1101-232, lES 0229H-200, and lES 0347-121 (see, 
e.g.. Fig. 8 in Fi nke et al. 2010). These unusual TeV spectral 
emission components are conventionally explained by (either 
leptonic or hadronic) emissions at the source, but they could 
also be explained by intergalactic cascade emissions. Non- 
simultaneous TeV excesses are also seen above th e extrapo- 
lation of the GeV fl ux in NGC 1275 dAbdo et al. 2009b) and 
the core of Cen A (lAbdo et al. 2010dl) . but because of their 
proximity, these excesses are unlikely to be UHECR-induced 
emissions made in intergalactic space. 

FigurelUdemonstrates that lES 0229H-200 can be fit by both 
the 7-ray induced cascade and proton-induced cascade emis- 
sions. Because of the uncertainty in EBL models, it is not easy 
to distinguish between the two possibilities at ^ 0.1 - 1 TeV 
energies. At higher energies, however, our calculations show 
that UHECR-induced cascade emission becomes harder than 
7-ray-induced cascade emission resulting from attenuation of 
hard 7-ray source photons, for a given EBL model. More 
importantly, the emission spectrum measured as a result of 
the injection of VHE/UHE photons at the source is strongly 
suppressed above ~' 10 TeV for a wide range of EBL mod- 
els, whereas a cosmic -ray-induced cascade displays a signif- 
icantly harder spectrum above this energy, and detection of 
> 25 TeV 7 rays from lES 0229H-200 is only compatible 
if the 7 rays are hadronic in origin. This is because UHE 
protons (and UHE nuclei) can inject high-energy pairs over 
the Bethe-Heitler energy loss length (Abh ^ (A/Z^) Gpc at 
Ea ^ A 10 eV) that is typically longer than the effective loss 
length of VHE/UHE photons. For steady, non-variable 7-ray 
sources, this intergalactic cascade signal induced by UHECRs 
provides a crucial probe of UHECR sources. Its identifica- 
tion would demonstrate that a distant blazar is an UHECR 
source through electromagnetic channels, which provides an- 
other important clue besides 7-ray variability. Identifying this 

' On the other hand, the emission of 7 rays and neutrinos produced in- 
side the source of primary UHE nuclei is hmited by the nuclear survival 
condition, as shown in Murase & Beacom (2010a; 2010b). Given that the 
observed UHECRs are dominated by heavy nuclei, this limitation is also 
appl ied to neutrinos produced outside the source, i.e., cosmogenic neutri- 
nos jMurase & Beacom 2010al) . 
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Fig. 8.— Spectral fits to HESS and VERITAS data of lES 0229+200. Blue 
data points are from HESS I Aharonian et al. 2007b). and red data points are 
preliminary VERITAS data (Perkins et al. 2010). The curves labeled "E20, 
low IR" and "El 9, low IR" are the cascade spectra initiated by the E^- in- 
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£max _ jq14 gY pf|Q(Qjjs \yj[[j ^ - 5^4 produced at the source for the low-IR 
EBL model. Double dot-dashed and dotted curves give, respectively, the Sc 
di fferential sensitivi ty for 5 and 50 hr observations with CTA (configuration 
E; IActis et al. 201 II) . 



feature by future Cherenkov detectors such as CTA or HAWC 
is poss ible, and the diff erential sensitivity goal of CTA is 
shown (lActis et al. 201 lb . Note that this is a differential sen- 
sitivity curve with the requirement of 5a significance for 50 
hour observations per bin, with 4 bins per decade. This is a 
much more stringent requirement than detection of a source 
with 5(7 based on integrated flux, which can be divided into 3 
data points with « 3ct significance each. Given the differential 
CTA sensitivity for a 50 hr observation, the spectral hardening 
associated with hadronic cascade development can be clearly 
detected. 

It is theoretically expected that cosmic -ray-induced and 
7-ray-induced cascade emissions are more easily discrim- 
inated in higher redshift sources. For the 7-ray-induced 
cascade, there should be a cutoff because of 77 pair cre- 
ation by the EBL, while spectra of the cosmic -ray-induced 
cascade emission are hardened by the continuous injection 
through the Bethe-Heitler process. Hence, deep observations 
at > TeV energies by CTA or HAWC for moderately high- 
redshift blazars will also be important to resolve this question, 
along with detailed theoretical calculations for individual TeV 
blazars. 

Now that IceCube has been completed, it has started to 
give important insights into the origin of UHECRs by itself 
and with GeV and VHE 7-ray observations. But detection 
of neutrino signals produced outside the source seems diffi- 
cult for high-peaked BL Lac objects and FR-1 galaxies, be- 
cause the point source flux sensitivity at > 10 PeV is order 
of- 10"" TeV cm-2 s"' (ISpiering 201 UlAbbasi et al720TTh . 
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9. — Spectra of UHE proton-induced neutrino emission for various 
redsliifts. The parameters used liere are the same as Figure 7. 



which is typically larger than the expected neutrino fluxes, as 
shown in Figure 9. On the other hand, the cumulative (dif- 
fuse and stacked) background neutrino flux may be detectable 
especiaUy for E™'' > 10^" eV (c.f. 'Anchordoqui elaLZOOl 
iTakami et al. 20091 and references therein), which is possi- 
ble in hadronic models with large magnetic fields in jets. 
For Zs"^" = 10' eV, however, protons mostly interact with 
the EBL, and the expected flux is lower than the proton 
case even if nuclei can be accelerated up to Ea = ZEp (e.g., 
[Anchordoqui et al. 2007). If BL Lac objects and FR-1 galax- 
ies are the main sources of UHECRs made mainly of ions with 
Ea < ZIO eV, the cumulative neutrino background would be 
difficult for IceCube to detect. 

Next, let us discuss the UHECR luminosity required to ex- 
plain such extreme TeV blazars in the intergalactic UHECR- 
induced cascade scenario. In Figure 8, with fi™^" = lO' eV 
and /? = 2, the inferred isotropic-equivalent UHECR luminosi- 
ties (at the source) are Luhecr — 10 ^ erg s"' (for the filament) 
and LuHECR — 5x10 erg s"' (for the cluster), respectively, 
while L^ ~ 10 erg s"' when primary 7 rays are injected. 

If no structured EGMFs are there, the required isotropic- 
equivalent UHECR luminosity (at the source) is Luhecr ^ 
10-10 erg s"' (i.e., the corresponding absolute, beaming- 
corrected cosmic -ray luminosity Luhecrj' ^ 10 ^ erg s"'), 
which is consistent with Essey et al. (2011), Razzaque et al. 
(2012), and the calculation for lES 0229H-200 shown here. 
It is also consistent with the proton power needed in hadronic 
models for typical, variable BL Lac objects. Where structured 
EGMFs play a role, the required UHECR luminosity becomes 
much larger In fact, we obtained Luhecr ^10 - 10 ^ erg s"' 
(for filaments) and Luhecr '^^ 10 ^ - 10 erg s"' (for clusters), 
depending on the EBL model and spectral indices, when we 



assume fi™" =10 
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eV. Such UHECR luminosities, obtained 
with the structured EGMFs, seem rather extreme, since the 
total cosmic-ray luminosity including low-energy cosmic rays 
is at least ^ 20 times larger for p >2 and fi™'" =10 GeV. In 
the cluster EGMF case, this means an absolute luminosity of 
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(13) 

where Mbh is the black hole mass. Therefore, the intergalac- 
tic UHECR-induced cascade interpretation becomes problem- 
atic, if runaway UHECRs are significantly isotropized and/or 
the spectral index of cosmic rays is steep enough. Such 
isotropization may be realized by some plasma instability, or 
the structured EGMFs and/or magnetic fields in radio bubbles 
or lobes accompanied by radio-loud AGN (see below). 

Note that our conclusion from Table 2 does not hold in 
the intergalactic hadronic cascade interpretation of extreme 
blazars since the SSC model is here abandoned. But one may 
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adopt £"" ~ 10 eV, motivated by results of the SSC mod- 
eling for typical, variable blazars (see Table 2). On the other 
hand, higher /i™" is also possible and a proton spectrum with 
higher E™" is favored in view of smaller deflections and re- 
laxed luminosity requirement to fit TeV data (see Figure 6). 
However, similarly to the proton synchrotron blazar model 
for variable BL Lac objects, the proton synchrotron compo- 
nent is expected as well as the intergalactic hadronic cascade 
component. 

It is useful to compare those luminosities with the required 
UHECR energy budget indicated from UHECR observations. 
From recent PAO observations, the local UHECR energy bud- 
get above lO'^'^ eV is a few x lO"*"* erg Mpc""' yr"' . For the lo- 
cal blazar density, n,, -- lO"*"^ Mpc"^ jPadovani & Urrv 1990h . 
the inferred isotropic-equivalent UHECR luminosity is 
^UHECR ^10 ^ erg s"' (regarding blazars as radio-loud AGN 
pointing toward us). This is much smaller than the cosmic- 
ray luminosity required for explaining extreme TeV blazars, 
and implies that those distant radio-loud AGN with hard VHE 
spectra should be rarer and more powerful in cosmic rays than 
nearby AGN responsible for the observed UHECRs. For lES 
0229H-200, the single-source flux is ^ 10 % of the observed 
UHECR flux so that the anisotropy can be used as a useful 
probe. 

We demonstrated the importance of structured EGMFs that 
help isotropize the trajectories of UHECRs, though the EGMF 
strengths are still uncertain. In addition, there are other causes 
that can diminish the beaming of UHECRs and resulting cas- 
cade fluxes. One arises from plasma instabilities induced 
by cosmic rays (Murase et al. 2011, in preparation). Sec- 
ond, radio lobes of powerful radio-loud AGN, like in the 
case of Cen A with B ^ 1 fxG, would also isotropize UHE- 
CRs jDermer et al. 20091) . as might radio bubbles from the 
jets of typical FR-I radio galaxies and aligned counterparts. 
These magnetic fields seem relevant in order that cosmic rays 
from relativistic jets of radio-loud AGN to contribute to the 
observed flux of UHECRs. Indeed, for nearby radio-loud 
AGN, the UHECRs must be significantly isotropized, since 
there is no blazar (i.e., aligned radio-loud AGN) within ^ 
100 Mpc'" and no evidence of cross-correlation with nearby 
blazars such as Mrk 501 and Mrk 421 (Dermeret al. 2009). 

'" In other words, the "apparent" UHECR source den- 
sity indicated from analyses of auto-coiTelation satisfies n^ > 
Ur^ Mpc"3 jKashti & Waxman 20081 : ITakami & Sato 20091) . which is 
larger than the local blazar number density, «i ~ 10"*' Mpc"'. If UHECRs 
are isotropized rather than beamed, one may compare it to the local FR-1 
galaxy density, «., ^ 10^ Mpc"' (e.g., Padovani & Urry 1990), which is con- 
sistent with the lower limit on the apparent UHECR source density. Then, the 



10 — 10 A/q black infened UHECR luminosity per source is typically LuHECR.; ~ 10'*' erg s 
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The isotropic-equivalent UHECR luminosity (at the source) 
^UHECR ^10 erg s"' at ^ 100 Mpc will lead to overproduc- 
tion of the observed UHECR spectral flux. 

4. DISCUSSION AND SUMMARY 

In this work we studied BL Lac objects and FR-I radio 
galaxies as potential UHECR sources in light of recent Fermi 
and imaging atmospheric Cherenkov telescope observations, 
and considered how future CTA, HAWC, and other high- 
energy 7-ray experiments might test the origin of the 7 rays 
from this class of blazars. 

If one accepts the standard synchrotron/SSC model for typ- 
ical, FR-1 galaxies and highly variable BL Lac objects that 
comprise the majority of VHE radio-loud AGN, the proton 
maximum energy is typically ^1-10 EeV unless UHE pro- 
tons are produced as rare transient events, and only heavier 
nuclei normally reach the > 10^" eV energies. In terms of the 
maximum energy, a heavy-ion dominated composition can be 
compatible with the standard SSC model because Fe nuclei 
can be accelerated to > 10^"^ eV while surviving against pho- 
todisintegration (if 5 > 20; see equation 8). An open issue of 
the heavy-ion dominated composition scenario of radio-loud 
AGN is how the significant amount of heavy nuclei is loaded 
in AG N jets, which i s suggested from the R40 composition 
results ("A braham et al. 2010i) and the observed isotropy in ar- 
rival distribution at ~ lO'^ eV Z{\Ea.2Q.5 (lAbreu et al. 201 Ih . 

On the other hand, if hadronic models are adopted for 
typical, FR-1 galaxies and highly variable BL Lac objects, 
then the observed VHE emission from these objects could 
be proton synchrotron radiation if protons are accelerated 
up to ^ 10""^ eV, which requires strong magnetic fields, 
B' ^ 10-100 G that could be found in the inner jets of the 
radio-loud AGN. Such hadronic models can be compatible 
with a proton-dominated composition. Especially for lumi- 
nous blazars with spectacular flares and low-peaked BL Lac 
objects with scattered radiation fields, one may expect high- 
energy neutrino s produced in inner jets a s one of the hadronic 
signatures (e.g. JAtoyan & Dermer 20011) . 

In either of the synchrotron/SSC or hadronic model, we 
mainly considered the blazar zone in the inner jet as the emis- 
sion region of 7 rays. However, some recent studies based on 
simultaneous radio and 7-ray observations are questioning the 
standard i dea that the blazar reg ion is located near the AGN 
core (e.g., iMarscher et al. 2010h . For example, in the case of 
3C 345, Schinzel et al. (2012) proposed that the emitting re- 
gion is located at ^ 23 pc along the jet. 

Extreme TeV blazars, in sources like lES 0229H-200, 
lES 0347-121, H 23 46-309, and lES 1101- 
232 (iNeronov & Vovk 201 Oh . are extreme both in their 
deabsorbed TeV spectra and their quiescent, non-blazar-like 
behavior. Their hard source spectra can be explained by 7 
rays that are produced either via electronic SSC or hadronic 
processes in inner jets, but hard EBL-deabsorbed VHE spec- 
tra typi cally require extreme so urce parameters or a specia l 
setup (iTavecchio et al. 201 lb IZachar opoulou et al. 201 Ih . 
An alternative interpretation for the extreme blazars, whose 
variability is slow or absent, comes from intergalactic cas- 
cade emissions. 7-ray-induced or UHECR-induced cascaded 
emissions can make slowly variable or almost non-variable 
components provided that the IGMF is weak enough, and 
the former may be seen as a slowly variable pair-echo 
component, following a more rapidly variable component. 

We examined these possibilities with numerical calcula- 
tions, taking into account effects of structured EGMFs in 



filaments and clusters, and demonstrated that the structured 
EGMFs would play an important role on the proton-induced 
cascade emission, and that cosmic rays are significantly 
isotropized especially for the maximum proton energy of 
< 10 eV In this case, rather large cosmic-ray luminosities 
are required in order to explain the emission from extreme 
TeV blazars such as lES 0229H-200 by the cascade radiation 
induced by UHECRs. Note that adopting £;f" - lO'"* eV 
is motivated by values obtained from the synchrotron/SSC 
modeling for typical, variable radio-loud AGN rather than ex- 
treme, non-variable blazars. One should keep in mind that the 
intergalactic cascade scenario itself does not tell about /i™". 

Hence, one may consider protons with /s™^" ^ 10"" eV, where 
strong deflections in the structured regions can be avoided. 
As a result, the required UHECR luminosity to power VHE 
emission from extreme TeV blazars becomes more reason- 
able, which is also comparable to values needed by the proton 
synchrotron model for typical, variable TeV blazars. 

There are some issues in the UHECR-induced cascade sce- 
nario that need more discussion. Even if cosmic rays are 
completely beamed upon entering the voids of intergalactic 
space, the distant and luminous extreme blazars should be 
more powerful cosmic-ray sources than the nearby weaker 
FR-1 galaxies found within the GZK radius that are assumed 
to be responsible for the highest energy UHECRs. To avoid 
requiring excessively large beaming-corrected cosmic -ray lu- 
minosities to make steady TeV radiation, it is therefore bet- 
ter if the UHECRs from extreme TeV blazars are beamed. 
On the other hand, if BL Lac objects and FR-1 galaxies are 
the main sources of UHECRs, then the escaping UHECRs 
must be significantly isotropized, because most radio galaxies 
within the GZK radius have misdirected jets. This isotropiza- 
tion could be caused by radio bubbles or the lobes of radio 
galaxies. The contrary behavior may be allowed if the distant 
extreme blazars from which steady emission is detected are 
in some ways special, e.g., their host galaxy exists in a region 
with a very weak EGMF 

Furthermore, note that in the case of rare transient activities 
with short duration fdur and longer quiescent periods f, 



qui 



> 



max[ArcR, AfiGv] between events, the intergalactic cascade 
scenario for extreme blazars may not work. The timescale 
during which cosmic-ray-induced cascade emission lasts for 
long time is increased by the factor max[ArcR, Aficvl/fdui- 
Then, the cosmic-ray-induced 7-ray flux is consequently re- 
duced by the cosmic-ray and pair-echo induced extension 
of this emission, and the required cosmic-ray luminosity 
has to be correspondingly increased, which could make ex- 
cessive demands on UHECR power. Such a rare tran- 
sient episode could, however, produce long-lasting though 
faint cosmic -ray-induced 7-ray emission that lacks associ- 
ated source emission, in analogous to TeV 7-ray-induced 
cascade radiation (e.g., iNeronov et al. 2010h . On the other 
hand, for repeating flaring activities with the time interval 
fqui < max[ArcR, AfiGv], the cosmic-ray-induced cascade 
emission can be regarded as almost persistent due to contri- 
butions from multiple flares. 

Observational tests of the properties of the extreme TeV 
blazars can reveal the radiation mechanism of BL Lac objects 
and FR-1 galaxies and provide a clue to UHECR acceleration, 
despite the above potential issues. First, continuing variabil- 
ity seai-ches with VERITAS, HESS, and MAGIC, and future 
studies with CTA are obviously important to determine if the 
emission must be made in the jet. In addition, discrimination 
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of the UHECR-induced intergalactic cascade from the 7-ray- 
induced cascade and attenuated source emission is possible 
from measurements at > 1 - 10 TeV energies. Detection of 
high-energy photons above 25 TeV from lES 0229+200 or 
above ^ TeV from more distant blazars, which may be real- 
ized by future 7-ray detectors such as CTA and HAWC, would 
be compelling evidence that this kind of object is a source 
of UHECRs. In addition to 7-ray observations, anisotropy 
searches with UHECR arrival directions provide another in- 
teresting test. 

In this work we focused on UHECR acceleration in the in- 
ner jets of BL Lac objects and FR-1 galaxies. Note that FSRQs 
and FR-11 galaxies are rarer but more powerful, and they can 
also be sources of UHECRs and neutrinos. Other scenarios 
such as shock acceleration at hot spots (e.g., Takahar a 19901 : 
[Rachen & Biermann 1993; Takami & Horiuchi 201 1 ) and co- 
coon shocks (e.g., iNorman et al. 1995t lOhira et al. 20101) 
are viable for these types of radio-loud AGN. But their 
local number density, n^ ^ 10"^ Mpc~ , appears to be 
too sr nall to avoid strong anisotropy in the local uni- 
verse (iTakami & Sato 20091) . In the heavy-ion dominated 
composition case, not only blazars a nd radio galaxies, but also 
radio-quiet AGN dPe'eretal. 2009h could be sources of UHE 
nuclei. 

In summary, we have considered observational implications 
of BL Lac objects and FR-1 galaxies as steady sources of 
UHECRs. Within the standard synchrotron SSC model for 
typical, variable BL Lac objects and misaligned counterparts, 
acceleration of UHE protons to energies > 10^** eV is unlikely. 



so the composition of higher-energy cosmic rays should be 
dominated by heavy ions within the framework of this model. 
The intergalactic cascade emission has to be sub-dominant 
for highly variable blazars and radio galaxies, while it can 
play a role on the spectrum of slowly variable or non-variable 
objects, especially extreme TeV blazars. If the TeV spec- 
trum of those blazars is produced by UHECR-induced cas- 
cade emission, then structured EGMFs, which can signifi- 
cantly isotropize protons, increase the luminosity demands on 
these sources. The intergalactic cascade emissions induced 
by VHE/UHE photons and UHECRs from a distant source 
can be distinguished by future multi-TeV observations from 
CTA and HAWC. In particular, detection of > 25 TeV photons 
from relatively low-redshift sources such as lES 0229+200 or 
> TeV photons from more distant sources would favor such 
objects as being sources of UHECRs. 
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